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Abstract. We report the detection of the H2 u=l-0 S(l) 
line at 2.122^m, from RNO 91 in the L43 dark cloud, 
which is known to be a T Tauri star surrounded by a 
1700AU disk structure (containing ices) and a weak out- 
flow. The non-detection of the H2 v=2-l S(l) line at 
2.247/im suggests shock excitation rather than fluores- 
cence. The emission is extended spatially up to 9" in the 
north-south direction. The line intensity peak (FWHM ~ 
3") corresponds to the star RNO 91 which is embedded in 
a cocoon of gas and dust. The observed H2 emission from 
this cocoon may be attributed to embedded Herbig-Haro 
like knots. The H2 line flux in the central 2" x 3" is esti- 
mated to be 7 x 1CP 14 ergs sec _1 cm~ 2 , which indicates a 
mass flow rate of 4 x 10~ 8 M Q yr _1 . Furthermore, narrow 
band image taken through H 2 1-0 S(l) filter is presented, 
which reveal a tilted disk and bipolar outflow structure 
that agrees with earlier observations and models. We show 
that this disk/outflow system is a unique case. 
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1. Introduction 

The near-IR molecular hydrogen emission lines are recog- 
nised as being important tools in studies of star formation 
(Shull & Beckwith 1989] ). The excitation of these lines in- 
volves mainly two competing processes: (i) shock heating 
and (ii) UV fluorescence (Burton 1992). ft is possible, how- 
ever, to distinguish between these two processes by mea- 
suring the ratios of intensities of lines arising from two 
different vibrational levels (Sternberg & Dalgarno 1989, 
Hora & Latter 1994). The origin of shocked molecular 
hydrogen emission from the spatially unresolved region 
close to the star can originate from outflows or from ac- 
cretion shocks in a disk, as in the case of Infrared Com- 



panions (IRC's) (Herbst et al. |l995j Koresko et al. |1997| ). 
RN0 91(a = 16 h 34 m 29.3 s ,<5 = -15°47'01")is one of the 
only two known PMS stars in the L43 dark cloud in 



Ophiuchus. It was classified as a MO. 5 type T Tauri star 
by Leverault ( 198SQ , based on an optical spectrum that 
showed strong Ha emission. An outflow driven by this 
star was identified at millimeter wavelengths and shown to 
have spatia lly se parated red-shi fted a nd blue-shifte d lobe s 
(Leverault |l988i Myers et al. |l987[ Bence et al. |l998|) . 



Ho wever , optical images and spectra obtained by Schild et 
al. ( 1989 ) showed that the outflow does not have any emis- 
sion indicative of shocked material around RNO 91. Opti- 
cal and infrared photo metry (U-band to L -band ) was ob- 
tained by Myers et al.( |l987| ). Heyer et al. (|1990[) obtained 
the JHK photometry and H-band polarimetry of this ob- 
ject which showed a disk type structure for the first time. 
Weintraub et al. ( 1994 , hereafter, W94) demonstrated by 
their K-band polarimetric image and 3-5 fim spectra that 
RNO 91 is surrounded by a disk-like structure of radius 
1700 AU comprising frozen H2O, CO and possibly XCN. 
They also showed that the polarization center does not co- 
incide with the intensity peak identified as RNO 91. Our 
speculation that this result indicates an IRC to RNO 91 
was quickly laid to rest by the Shift and Add imaging by 
Aspin et al. ( 1997] ) (hereafter A97), which does not show 
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any secondary source within a 3.7" square region. In this 
paper we present near-infrared spectra and narrow band 
images of the object in the K band region. We also report 
the detection of the H 2 v— 1-0 S(l) line at the source and 
discuss its implications. 



2. Observations and Data Reduction 

Near infrared spectroscopic observations were made on 
March 25, 1998 at Gurushikhar 1.2 m Infrared Telescope 
(GIRT), Mt. Abu, India. A Near Infrared Camera / Grat- 
ing Spectrometer based on a HgCdTe 256 x 256 focal plane 
array was used to obtain the observations. The grating 
spectrometer was used in a configuration that yielded a re- 
solving power of X/5X = 1000 with a l"/pixel plate scale. 
The slit was two pixels wide and oriented along the N-S 
axis. The atmospheric seeing and image motions were be- 
low 2" during the observations, which was measured from 
imaging data obtained just before the spectroscopic obser- 
vations. Data acquisition and reduction were done using 
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Fig. 1. Contours of H 2 2.122/mi line from RNO 91 spec- 
trum showing the extended emission. Notice the different 
lengths of emission in the N-S from the center that in- 
dicates the tilt of the outflow axis. Continuum is partly 
subtracted using an adjacent strip of the same spectrum, 
which has caused the artifact of width variability. 

standard procedures. We used the RNO 91 K band pho- 
tometric fluxes (0.48 Jy) given by Myers et al. ( 1987 ) for 
flux calibrating the spectrum. 

Narrow band (AA = 0.02/im) images through H2 
(A = 2.122/im), Br7(A = 2.165/im) and continuum (A = 
2.104/im) filters were obtained by the United Kingdom 
Infrared Telescope (UKIRT) Service Observing Program 
on September 8, 1998 using the facility near-IR imager 
IRCAM3. IRCAM3 employs a 256 x 256 InSb array; the 
optics used gives a pixel scale of 0.280". There was a small 
defocusing problem that occurred during these observa- 
tions due to variable seeing, resulting in image elongation 
in the N-E, S-W direction. These uncertainties were esti- 
mated to be about 1.2". Continuum subtraction was not 
carried out because of the focusing problems in these im- 
ages. It should also be noted that there is some ghosting 
that had occurred in these images. These are identified as 
ghost images since they occur at exactly the same posi- 
tion with respect to the main source both in object and 
standard star frames. 



3. Results 

The contour map of the spectrum shown in figure 1 dis- 
plays the extended H2 emission along the N-S direction 
(the slit axis) corresponding to the outflow axis. The in- 
tensity peaks at -3" and +2" represent H2 knots in the 
outflow. From gaussian fits to the continuum emission pro- 
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Fig. 2. Spectrum of RNO 91 integrated over 9 rows rep- 
resenting an area on-source of 2"x 9". 



file at different positions along the dispersion axis, we find 
that the FWHM of the continuum strip (measured N-S) is 
about 3.0", although the seeing on the night of observation 
was between 1.5-2.0". The extra width in the stellar con- 
tinuum strip could be attributed to scattered light from 
the cocoon surrounding the star. We have therefore ex- 
tracted the source spectrum by integrating the emission 
along three and then five rows, representing an on-source 
area of 2" x 3" and 2" x 5". The relative intensities of 
the photospheric Nal and Cal features remained the same 
in these two extractions. However, in the 3"-wide extrac- 
tion, the 2.122/im line intensity was considerably smaller 
relative to the photospheric features. This result confirms 
that the broadening of the continuum strip is indeed due 
to scattered star light. 

Figure 2 shows the spectrum of RNO 91 in the wave- 
length region 2.05/xm to 2.3/iin. The spectrum is obtained 
by integrating nine rows, covering 9" along the N-S slit 
axis centered around the star. The spectrum displays 
prominently the H2 v—1-0 S(l)line at 2.122/im and the 
photospheric Nal and Cal absorption features at 2.20/im 
and 2.26/mi respectively. The 2-1 S(l) line at 2.247/im is 
below the noise level. The spectrum also displays the Br7 
emission line at 2.167/im. These features are marked in 
the figure. Most of the other feat ures s een in absorption 
are telluric in nature (Chclli et al. 1997 ). 

The excitation mechanism for H2 emission can be in- 
ferred from the 2-1 S(l)/l-0 S(l) flux ratio (e.g. Luh- 
1998| ) . An estimated upper limit to the 2-1 
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S(l) line flux yields this ratio to be 0.16 or less, point- 
ing strongly towards shock excitation. However, flourcs- 



Nanda Kumar et al.: Shocked H2 emission from RN0 91 



3 




outflow lobe in the south can be seen clearly. These struc- 
tures are more evident in this narrow band image, than 



in the earlier K' image of Hoddap (1994) or the polari- 
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metric image of W94. W94 had shown the existence of 
frozen H2O, CO and possibly XCN in a disk structure 
of radius 1700AU around RN0 91. The existence of these 
ices on grains within this disk structure was proven by 
the absorption features found in a 3 — 5/Ltm spectra. Their 
K-band polarimetric map also conforms wit h mo dels of 
scattering from disks (Whitney & Hartmann 1992 ). W94 
suggest a scenario with a flared disk viewed at an angle of 
30° with the northern outflow lobe tilted away from the 
observer. The narrow band image shown in figure 3 (note 
specifically the patch of continuum emission to the south 
of RN0 91) clearly support this scenario. In addition, it 
can be seen from fig 1 that the H2 emission is extended 
more in the southern direction than in the northern di- 
rection which also confirms the tilt of the outflow axis, 
the northern lobe of the outflow being obscured near the 
source by the disk structure. 



Fig. 3. Contour plot of RN091 taken through narrowband 
filter, illustrating an east-west disk associated with the 
source, as well as nebulous emission to the south. The 
contours measure 5,10,15,20,40,80,160,320 and 640x the 
standard deviation to the mean background level in the 
image. 



cent excitation in a high density regime may also produce 
a "shock-like" 2-1 S(l)/l-0 S(l) ratio, because of ther- 
malisation of low-energy vibrational levels (Sternberg & 
Dalgarno 1989] , Burton et al. 1990p . In a high density re- 
gion one would still expect to see emission from the u=3 
level at 10% of the 1-0 S(l) line (rather than 1%, as is ex- 



pected in a shock; see e.g. Luhman et al. 1995). However, 



in our spectrum the 3-2 S(3) line is unresolved from Nal 
absorption, and the 3-2 S(l) line is outside the wavelength 
range. 

We resolve this issue by estimating the two most cru- 
cial parameters that decide the efficiency of the UV fluo- 
rescence, namely, the gas number density and the UV flu x 
scaling parameter \ ( see Sternberg and Dalgarno |1989D . 
We estimate an upper limit for the gas number density to 
be 6 x 10 5 cm -3 , considering an A„ ^9 (Myers et al. 1987 ) , 
an outer disk radius of 1700AU (W94), and an inner disk 
radius of 0.01AU corresponding to the dust evaporating 
radius for an M0 star. Note that this density is for a gas 
disk of 1700AU and the actual regions from where we ex- 
pect the H2 emission are well below a radius of 200-300AU. 
On the other hand, the UV flux scaling parameter \ f° r 
an M0 type star is much less than 1. From these, the UV 
fluorescence is expected to be of little significance for the 
excitation of H2 lines. 

Figure 3 shows a narrow band image through H2 filter, 
without continuum subtraction. The disk structure and an 



4. Discussion 

From the spectrum we measure an H2 flux of 7x 10" 14 ergs 
sec -2 cm -2 (integrated over 3 rows, representing an area 
on-source of 2" x 3"). If we attribute this flux to shocks, 
we can estimate the mass flow rate M using the relation 



eL = GMM/R = l/2Mv 2 



(1) 



where e is the ratio of the total energy in the shock to 
the strength of the 1-0 S(l) line, L is the observed H2 line 
luminosity, and v the shock velocity. We adopt a value for 
v of 30 km s , a value that is optimum for produc ing H 2 
line emission. We also assume that e = 50 (Smith 1995 ). 
Together, these yield an estimate for mass flux of M = 
4 x 10" 8 M yr" 1 . 

As demonstrated in the previous section, the relative 
change in the intensity of H2 from that of the photospheric 
features in the spectra extracted with different widths of 
the continuum strip shows that the H2 emission close to 
the star likely originates from a region different to that 
traced by the photospheric lines. The fact that the H2 
line emission is extended along our N-S slit strongly sug- 
gests that the H2 is associated with an outflow. By ob- 
serving line emission coincident with the RNO 91 stellar 
continuum, we may therefore be tracing the outflow all 
the way back to the source. The total flux measured from 
the observed extended emission (integrated over 9 rows 
representing an area on-source of 2" x 9") is 1.5xl0" 13 
ergs sec" 2 cm" 2 . This represents an average flux dis- 
tribution of 3.5 xlO" 4 ergs sec" 2 cm" 2 sr~ 2 . Assuming 



Lh2 ~ 10 x L S n\ (Davis & Eisloffel 1995), the extinc 



tion corrected S(l) line flux yields a net H2 luminosity of 
0.002L©. Using the estimated value of M and the kine- 
matic age of the outflow, the H2 luminosity can be shown 
to represent a net warm H2 mass of ~ 5 x 10~ 4 M Q , where, 
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a 9" long flow with a velocity of 30 km s _1 is used to cal- 
culate the age of the outflow. 

The H2 line flux measured from within the central 3" 
could be associated with the near-IR counterparts of HH 
knots embedded within the cocoon that are excited by the 
outflow and accretion shocks in a disk. However, excitation 
of all of the observed "on-source" H2 emission due to infall 
is not feasible, since the mass accretion rates derived from 
such infall will produce a K-band extinction (Herbst et 
al. 1995) that is an order of magnitude higher than the 
measured value by Myers et al. ( |1987| ). But some fraction 
of the observed emission could originate from accretion 
shocks, with the rest from the outflow, as in the case of 
T Tauri. It is interesting to note that the high resolution 
images obtained via Shift and Add imaging at UKIRT by 
A97 reveal a nebulous feature which is about 1" north 
of the star. This could be a shock excited feature in the 
outflow. Alternatively, this feature - in conjunction with 
the fainter features to the S-W of the star - appears to 
be part of an ellipse whose major axis is ~ 200AU as 
measured by us using the published images of A97. In 
fact these features may be visualized in terms of a disk 
tilted by 30° to the north. The absence of any feature to 
the south may be attributed to the obscuration of the disk 
by the outflow. These arguments support the W94 data 
and their model, as do our narrow band images shown in 
Figure 3. 

While gas and dust disks are known to exist around 
PMS stars (Beckwith et al. |l990| , Dutrey et al. |l994j ), there 
exists a small sample of objects with disk structures that 
contain icy mantles where ices are found on dust particles 
in the protostellar envelopes (see Chiar et al. 
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and ref- 
erences therein). RN0 91 possesses by far the largest icy 
gas/dust disk structure (20") and also shows shocked H2 
originating in an outflow and possibly an accretion disk. 
These features make RNO 91 unique and an interesting 
object for further studies in millimeter and infrared wave- 
lengths. 

5. Conclusions 

An infrared spectrum of RNO 91 shows emission of 
shocked molecular hydrogen from an outflow in the N- 
S direction and from a spatially unresolved region close to 
RNO 91. We estimate a mass flow rate ofAf = 4xl0~ 8 
M© yr _1 , based on the line fluxes from this spatially un- 
resolved region around RNO 91. The line fluxes from the 
spatially extended outflow yield Lh2 ~ 0.002L Q repre- 
senting net warm H2 mass of ~ 5 x 10 -4 M©. The outflow 
seen here in H2 emission, extending roughly N-S, appears 
to support the tilted disk + outflow model of W94, where 
the northern flow lobe is tilted away from us at an angle 
of 30° to the plane of the sky. Our narrow band images 
also support this scenario. We suggest that the H 2 flux 
from the spatially unresolved region around the source 
could originate in "HH-type" knots embedded in the co- 



coon surrounding the star. However, accretion shocks can- 
not entirely be ruled out based on our observations. We 
argue that RNO 91 is a unique disk/outflow system. 
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